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1. Introduction 

Abstraction is a common tool in physics, serving the better comprehension of complex 
phenomena. The most evident example is the atom, encompassing an intricate 
underlying structure of electrons and hadrons, the latter themselves composed of quarks 
and gluons. This internal structure can of course be completely neglected while 
discussing the properties of gases in statistical physics. Such multi-level abstraction 
is especially common for solid-state physics, where excitons are formed from electrons 
and holes, themselves being complicated elementary excitations formed from several 
electronic levels of atoms constituting the solid. The excitons can be described as 
massive quantum particles, which can couple with photons and form new particles of 
an even higher level of abstraction - composite bosons called exciton-polaritons. 

When a Bose condensate [1] is formed, one can consider its weak excitations as 
elementary particles (called Bogolons), forgetting the nature of the underlying bosons. 
But the weak elementary excitations are not the only type of interesting perturbations 
which can occur in a Bose condensate. The topological defects |2] resulting from the 
non-linearity of the Gross-Pitaevskii equation and the irrotationality of the macroscopic 
condensate wavefunction are currently (and since quite a long time) in the focus of 
intense theoretical and experimental research and will be the main topic of the present 
work. We shall see how their behavior can be described in terms of relativistic "material 
points" and "point charges" (an easy way) or in terms of underlying local spin dynamics 
(a harder way). The importance of such analogies as "magnetic charges" will thus 
become especially clear. 

In a recent paper [3], we have discussed the behavior of single half-solitons|l] and 
half- vortices [6] in spinor polariton condensates in the presence of an effective magnetic 
field (induced by the energy splitting between polarizations in wire and planar cavities). 
We have shown that such topological defects behave as "magnetic charges", accelerating 
along the applied magnetic field: that is, like analogues of magnetic monopoles. Such 
a behavior is a special feature of half-integer topological defects, which also applies to 
the so-called oblique half-solitons [7]. We have shown as well, that their stability is 
maintained thanks to the spin anisotropy[9] of polariton-polariton interactions. 

In the present work we analyze mode in details the behavior of integer topological 
defects in the presence of an in-plane effective magnetic field. We show that they 
become unstable and split into half-integer ones, which will separate and accelerate in 
real space due to the interaction with the field. We concentrate on the practical ways of 
their experimental excitation in wire and planar microcavities in a controlled manner, 
in order to study their eventual separation and we discuss the combined extrinsic and 
intrinsic mechanisms that favor the separation. 
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2. Effective magnetic fields in the Gross-Pitaevskii equation 

While the topological defects in Bose-Einstein condensates (single or multiconiponent) 
are relatively well-known, their spin texture and their interaction with magnetic field 
have been much less studied (see ref.|5]). First of all, let us see how an intrinsic effective 
magnetic field appears in the spinor Gross-Pitaevskii equation: 

^^-Qf = ~2^^^+ + ailV'+lV+ + a2|V'-lV+ (1) 

ih-^ = -i^M- + «il^-lV- + «2|V^+|V- (2) 

Here a\ is the interaction constant in the triplet configuration (parallel spins), and 
02 is the interaction constant in the singlet configuration (opposite spins), usually 
negative [11] and much weaker than the former due to the fact that it is a second- 
order process involving dark exciton states. These constants are usually calculated 
using the formula p^ a\ = 6xEi,a'^/S, where S is the normalization surface, a^ is 
the exciton Bohr radius, x is the excitonic fraction of the polariton state, and E^ is the 
exciton binding energy, while 0^2 is usually assumed to be 02 = —(0.01.. 0.2)^2, or simply 
neglected. Below, we are going to check the influence of this choice on the behavior of 
the condensates. 

The spinor Gross-Pitaevskii equation can be rewritten in the following form: 

+ ^^(i^.r-i^-nv^. (3) 

■k^^- ^^ A / , ^1+^2 /, , ,2 , I , |2X , 

-^^(i^.r-i^-n^- (4) 

It is easily noticed that the last term enters the two equations with opposite signs, 
and therefore can be described by the Pauli matrix a^ in the Hamiltonian. It can 
therefore be considered as an effective magnetic field along the z direction. This field is 
responsible for the so-called self-induced Larmor precession of the polariton pseudospin 
|12j . It is very important to note here, that the direction of this field is opposite to 
the ^-projection of the polariton pseudospin. Therefore, if the condensate is completely 
circularly polarized with ■?/'+ = y/n, ip_ = 0, n being the total density, the effective 
magnetic field will point in the negative direction of the z-axis: 

a = -^^(l*+r-W.-r) (5) 

In the polaritons case, where 02 is around ten times weaker than ai but negative, the 
effective field is strengthened and tends to lock circularly polarized states (strong density 
imbalance) in the system. It provides the natural stability of half-integer topological 
defects that carry particles with spin up (o"+) or down ((T_) at their core. In usual spin 
isotropic atomic condensates ai ~ 02, and therefore Q^ — 0- 
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Another effect, which is well known in planar or wire microcavities, is the splitting 
between the orthogonal linear polarizations at A; = (as opposed to the wavevector- 
dependent TE-TM splitting). This splitting is especially well resolved above the 
condensation threshold [13j thanks to the sharpening of the emission lines and induces 
an in-plane effective magnetic field pointing in a well defined direction, determined by 
the orientation of the cristallographic axes of a planar cavity or the orientation of the 
microwire|14j. Its interpretation in terms of an constant effective in-plane magnetic field 
is even more direct than in the case of the interaction-induced field in the z direction: if 
the X and y polarizations have different energies when the Gross-Pitaevskii equation is 
written on the xy basis, this splitting transforms into a term —hQiT'^P^/'^ in the circular 
polarization basis with the usual coordinate transformation rules tp± = (ip^ T ii^y) /v^- 

For a homogeneous wavefunction, the kinetic energy term is absent, and one 
can easily see that the spinor Gross-Pitaevskii equation for the wavefunction can be 
converted into an equation for the dynamics of the pseudospin vector S with the 
components given by: 

Sy = ^ (^-^;) (6) 

S.= (|^+r-|^-r)/2 

The equation for the pseudospin dynamics is very simple. It is obtained from the 
equation for the wavefunction by multiplying the corresponding equations by the 
wavefunction or its conjugate and taking the real or imaginary parts, in order to obtain 
the X and y projections. In fact, it is the usual equation for the precession of a magnetic 
momentum vector in a magnetic field: 

%^&.S (7) 

However, the field itself depends on the magnetic momentum via the non-linear term 
Cl = VLltUx ^ '^ S^U;, (8) 

3. Spin-topological defects in a magnetic field 

The scalar ID Gross-Pitaevskii equation is known to possess grey soliton solutions [1] 

(9) 



ijj [x — vt) = A/n 




iV2 



vt I v^ 
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here v is the speed of the sohton, ^ = h/^/2m,an is the heahng length, and c = \l anjva 
is the speed of sound. Let us consider the time evolution of the simplest extension of 
this scalar solution to a vectorial condensate: a half-soliton. In this part of the text, we 
will neglect the small interaction between opposite spins a^ for the sake of simplicity. 
The dark half-soliton solution for which f = reads [7]: 

V<, = y|tauh(^).V.- = y| (10) 
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where we have assumed no phase difference between the two components. This solution 
is shown in the figure 1(a) together with its pseudospin pattern. One can note that 
the half-sohton is a domain wall between the two regions of perpendicular linear 
polarizations (opposite pseudospins). We remind that the in-plane pseudospin field 
^ii = {Sx, Sy)'^ describing linear polarized states, makes a double angle with respect to 
the polarization angle. We see that the pseudospin pattern is divergent in that case, 
resembling the field of a point (magnetic) charge. The sign of the charge is given by the 
pseudospin pattern of the half-soliton (which is not necessarily divergent) and is defined 
by the relative phase between the two components. Let us first try to understand 
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Figure 1. (a) The two circular components cr+ (blue line) and a- (purple line) of the 
half-soliton wavefunction with the corresponding pseudospin vector field (red arrows), 
(b) Precession of the pseudospin under the action of the effective field H.. 



from "microscopic" considerations what will happen to this object when an in-plane 
effective magnetic field Qlt = ^x^x is applied. For this, one can consider the pseudospin 
dynamics at each point, since the two other terms in the Gross-Pitaevskii equation (the 
kinetic energy and the interaction energy described by the blue shift) compensate each 
other at t = 0. For the regions far from the soliton core |x| ^ ^, the pseudospin is aligned 
or anti-aligned with the magnetic field (which contains only a x component) and there is 
no evolution: dS/dt = 0. For the core region, there is a non-zero pseudospin projection 
on the z axis, and therefore the pseudospin will rotate (precess) around the magnetic 
field (which now contains both x and z components). Let us consider the initial moments 
of this rotation for the pseudospin in the center of the soliton: S{x = 0,t = 0) = Uz, 
while the total magnetic field is 



fl{x = 0, t = 0) = QxUx — 



ai — 0^2 
2h 






:iii 



Rotating around the (positive) fix component, the pseudospin (initially negative along z) 
gains a positive 5*^ projection and starts to precess around the positive flz field, turning 
towards the negative direction of the x-axis [see Fig. 1(b)]. This is the main result of our 
qualitative vectorial consideration of the polarization dynamics: the pseudospin in the 
center of the soliton gains a negative x-projection. Therefore, the domain of negative x 
pseudospin projection becomes larger, the domain of positive x-projection smaller, and 
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the wall between these domains is moving to the right. However, once the soliton core 
starts to propagate, the kinetic and interaction energy terms are no more compensated 
everywhere, and one cannot discuss the evolution of the system using the qualitative 
arguments based on polarization dynamics. At first sight, one may even think that the 
nonlinear system in question can be solved only numerically. 

However, an important insight into the behavior of the system can be gained by 
"changing the zoom". Forgetting about the complicated internal structure, a vectorial 
grey soliton in a Bose-Einstein condensate can be considered as a particle |H] with a 
negative effective mass (at least at low velocities). Moreover, the pseudospin pattern of 
this particle is the same as the field of a point magnetic charge in ID. The magnetic 
energy of the system can be found from the Hamiltonian as the usual scalar product of 
the field and the spin, and this magnetic energy depends on the position of the soliton 
because of the finite system size. Thus, one can evaluate the force acting on the magnetic 
charge from the magnetic field as a gradient of the magnetic energy with respect to the 
position of the soliton. This force will consequently accelerate the soliton. 

Considering the soliton as an elementary particle without internal structure means 
passing to the limit L ^ ^, where L is the system size (for example, the length of the 
wire cavity, which is usually of the order of 100/im). The healing length of a polar iton 
condensate for a reasonable blue shift of 1 meV expected for GaAs or CdTe cavities 
and a polariton mass of 5 x 10"^ of a free electron mass is ^ = 0.8/im. In this limit, 
the wavefunction of the soliton at t = becomes simply -ip^ = ^^/n/2siga.{x — xq), 
ip- = yn/2, where xq is the position of the soliton. In general, the tank function is 
replaced by the sign function. The magnetic energy of a condensate containing a HS in 
an external in-plane magnetic field is 

E^ag = - fi^/2-S)dx (12) 

Here S{x — Xq) = nsign(x — Xo)ux {xq is the soliton position) in the limit we consider, 
giving 

Eraag {x - Xq) = HQ^U (x - a;o) (13) 

dEma 
^ mag 3 Tl\Lx V-'-^J 

The force in Eq.(14) is therefore acting opposite to the direction of the magnetic field, 
but the acceleration will occur in the direction opposite to the force, because the mass 
of the soliton (which is a density notch in the condensate) is negative, at least at low 
velocities. 

For a grey soliton propagating at speed f , the phase shift induced by the soliton 
in the aj^ component is /S.0 = 2arccos(f /c) < n and the pseudospin projection Sx is 
reduced, which can be expressed as a renormalization of the magnetic charge. The 
correction to the charge is found as g = go(l — v^/c'^) (go is the charge at rest for a dark 
half-soliton) by integrating the solution ([9]) in the limit L ^ ^. The total correction for 
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the mass of the sohton and its charge gives the equation of motion 



uQt / v"^ 



2\ 3/2 



a = qo—^[l--] (15) 

mo V c J 

the same as in relativistic physics, integrating which one obtains 

t;(t) = ctanh('^^ty (16) 

assuming zero initial velocity. This trajectory is perfectly confirmed by numeric 
simulations as shown in Ref.p]. 

Thus, using the abstraction of a point magnetic charge proves particularly useful, 
because it allows solving the nonlinear spinor Gross-Pitaevskii equation analytically and 
gives a good qualitative understanding of observed phenomena. In the following we will 
discuss the soliton interaction and how to realistically generate currents of half-solitons 
and half- vortices within semiconductor microcavities toward "polariton magnetricity" , 
similarly to the "magnetricity" reported in spin ices [T7], but with much higher flow 
velocities (currents). 

3.1. Soliton and half-solitons interactions 

Dark (grey) solitons are the solution of the ID Gross-Pitaevskii equation provided that 
the condensate is formed of particles that repel each other («i > 0). Interactions 
between solitons themselves have been thoroughly analyzed an it is now well known 
that dark solitons tend to repel each other as well [15]. Such a behavior is illustrated 
in the Figj2] showing the degree of circular polarization and where we have solved the 
scalar Gross-Pitaevskii equation with the initial condition given by two dark solitons 
spatially shifted by 0.5 /im [see Eq.( l20l) below]. The solitons repel at the initial moment 
and then demonstrate a linear trajectory, which is a clear signature of their short range 
interaction. 

Now, what about the half-soliton in two-component condensates? Let us remind 
that a half-soliton corresponds to a soliton occurring in only one of the two components. 
Additionally, if the intercomponent interaction is not neglected, the presence of a half- 
soliton in one component obviously perturbs the other one in a fashion depending on the 
type of interactions, given by the sign of 0^2. Indeed, for q;2 > (0^2 < 0), the presence 
of a density minimum in the component carrying the soliton is seen as a potential well 
(barrier) by the other one and consequently a density maximum (minimum) appears at 
the position of the soliton (this effect is discussed in details with an analytical estimate 
of the density perturbation in Ref. |7J). Putting two half-solitons lying in different 
components next to each other induces a short range interaction within each component, 
between a soliton and the density maximum (minimum) induced indirectly by the other 
soliton. The half-soliton having a negative effective mass are repelled by the potential 
wells they encounter (which, from their point of view, are barriers). This argumentation 
leads to the following statement: For a^ > Q (02 < 0), half-soliton attract (repel) each 
other while they obviously don't interact for 02 = 0. This trend is shown in Fig. 2(b) and 
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Figure 2. Short range interaction between solitons initially separated by 0.5 fiva. 
(a) Integer soliton repulsion the colormap shows the density n = jV'l of the single 
component of a scalar condensate, (b) Half-solitons attraction, which shows dipolar 
oscillations for a2 — +0.2ai (c) half-solitons repulsion for a2 — — 0.2ai, the colormap 
in (b) and (c) shows the degree of circular polarization of the spinor condensate: 
Pc == ("-+ - n^)/{n+ + n_). 



(c). We note, that in the repulsive case, dipolar oscillations occur between the solitons, 
this is due to the effect of their mutual attraction and the slight initial displacement 
(the system can be described by a simple pendulum in this case). In the framework of 
polariton condensates where 02 < 0, the repulsion of half-solitons is therefore a strong 
asset for the separation of integer soliton into a pair of half-solitons of opposite charges. 
We note that bright solitons have also been described for exciton-polaritons [TU [19] , 
but these are not in the focus of the present work, although the same considerations of 
the interactions with magnetic field would apply to them as well. 



3.2. Solitons imprinting using Gauss- Laguerre beams 

In a recent work [20] we have made a proposal for the creation of half-integer (valid 
also for integer) solitons in quasi- ID polariton condensates under continuous quasi- 
resonant pumping. It consists in applying pulsed potential either equally to both circular 
components (integer topological defects) or to a single one (half- integer defects). This 
could be experimentally realized by sending short optical pulses at detuned frequency 
with respect to the main pumping laser. Such a setup would however require the use 
of two lasers and the cw-resonant injection would tend to imprint the phase of the 
condensate at any time under the pump spot, which would be harmful to the solitons 
characterized by local phase dislocations. Another method of hydrodynamic nucleation 
of solitons has been recently proposed theoretically in Ref. [2T] and then implemented 
experimentally [22] . However, this method requires the use of a defect, whose properties 
cannot be controlled at will. Here we propose a simpler configuration for the creation 
of half and integer solitons in quasi-lD microwires, while the 2D case will be discussed 
in the next section. 

Vortices have already been created in polaritons condensates by several means, 
among which is found the artificial phase imprinting induced by a probe carrying an 
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angular momentum: a " Gauss-Laguerre (GL) beam" [23J. Such a beam is obtained from 
a usual Gaussian beam scattered on a hologram containing a forklike dislocation. In 
numerical simulations, this type of pumping can be described by the following complex 
function 

Pgl (r, t, I) = AcWi^-^^y + iy-^yy 

fx-Ax\'2 ( y-^v Y ftzAlY ;. ■ . 

X e-\^^) e V -« ; e~y~^> e'^'^e-*^^^*, (17) 

written here in Cartesian coordinates. The exponentials give the " Gauss" part of the 
function, while the square root is the 1'^* order Laguerre polynomial, giving zero density 
in the center of the beam. I is the integer winding number that is to be transferred 
to the vortex state. The laser frequency ujgl should be slightly blue-detuned from the 
bare polariton mode, in order to make use of the bistability effect and obtain an almost 
fiat density profile of the condensate except in the center, due to the saturation of the 
pumping efficiency on the upper bistability branch [24]. 

As one can see from the sixth factor of Eq. flT7|) . where is the polar angle, the 
phase is changing continuously from to 2/7r encircling the beam center. In a pure ID 
system, the notion of angular momentum vanishes and grey solitons embody elementary 
topological excitations in BECs in place of vortices. The action of a GL beam is therefore 
to induce a local phase step of — /vr or smaller, equal to the one resulting from a cut 
of a vortex by a plane. It is illustrated in Fig|3]^a) showing a density slice of the GL 
beam together with its phase. This technique constitutes an efficient mean for soliton 
engineering in ID condensates. We note that, in order to allow the soliton to evolve 
freely, it is necessary to use a pulsed GL beam, which will also form the background 
condensate for sufficient pump intensity. 

Besides, experimental creation of half-solitons with independent selection of their 
phase (given by /) would require to separate a linearly polarized input laser in its two 
circularly polarized components (cr-t) using polarizers, to make one component (for a 
single half-soliton) or both of them (for a pair of half solitons) scatter on a hologram, 
and to recombine them on the sample. For the case where both components carry 
an angular momentum it looks difficult to recover a perfect spatial overlap of the two 
beams, which will naturally trigger the separation of the half-solitons, especially for the 
normal case of q;2 < 0. Additionally a relative phase 0o can be introduced between the 
two beams increasing the optical path of one component with respect to the other, ^q 
has a crucial impact on the pseudospin textures of HSs and therefore on their interaction 
with effective fields as follows from Eqs.(7) and (12). A scheme of an experimental setup 
is shown in the Figl3]^b). 

If we consider now a more realistic system such as a microwire [see FigjS] , we should 
take into account the transverse width of the sample and treat the case of a quasi- ID 
condensate for which the angular momentum cannot be neglected anymore. The central 
position of the pump spot crucially impacts the symmetry of the imprinted flow. We 
have performed numerical simulations of a 2 /im wide and 100 /im long wire using 
polariton parameters to highlight this feature. On one hand, if the spot is transversally 
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centered, a static and therefore dark soliton is nucleated in the wire, almost as in the 
pure ID case [see Figj3|^c),(d)]. On the other hand, a shift of the GL beam along the y- 
axis induces an extrinsic uniform propagation of the resulting grey soliton in a direction 
depending on the transverse shift and on the sign of the imprinted winding number I 
as it is illustrated in the figure [3]^c)-(f). A positive (negative) y-shift induces a soliton 
propagation to the right (left) for / < and reciprocally for / > 0. This effect can 
be qualitatively seen as a "rolling" of the particles on the boundaries, introduced by 
a gradient of angular velocity in the transverse direction. Moreover, the closer is the 
spot center to a boundary, the faster the soliton moves along the wire. This "extrinsic" 
motion has to be taken into account for an experiment, in which imposing a perfectly 
centered spot is elusive. Even in the simulation [see FigJSt^d)], where the precision on 
the centering is about 0.1 ^m, we observe a very small drift of the solitons. 

For the simulations, we use the typical parameters of modern GaAs cavities, 
including their outstanding lifetimes of about 30 ps [16]. We solve the spinor Gross- 
Pitaevskii equations for the photonic 0(f, t) and excitonic fields x(^) t) (coupled via 
the Rabi splitting Vr = 15 meV), fully taking into account the polarization a± = ±, 
the finite lifetimes r^ = 25 ps and t^ = 300 ps and the injection of the particles via 
^gl(^'^'^±''/'o): 

h"^ . . Vr ih , „, ^^ , ho.^ 



dt 2mA, ^^ ' 2 2ta, 



^^^ = - 7^—^4>± + -^x± - — 0± + Pgl + U(^± - 



■h^^^ ^^ A , ^^A *^ , ^ I |2 , I |2\ na\ 

^^-Qf = - 2^ '^^ T*^^ ~ 2^^^ ("i|x±l + «2|XtI j X± (19) 

X. X. 

Here m^ = 3.6 x 10~^r7io5 ""^x ~ 0.4?7io and mo are the cavity photon, the quantum 
well exciton and the free electron masses respectively. U{r) is a potential (e.g. the 
confinement potential for quasi-lD wires, or a wedge potential, see next section). 

3. 3. Natural acceleration of half-solitons in microwires 

Armed with efficient means of creating half solitons, let us now argue on their potential 
propagation and acceleration. We shall consider the following situation: A pair of 
half-solitons is created by the separated GL beams. They are slightly spatially shifted 
either initially or naturally due to an intrinsic noise (checked) that breaks the symmetry 
between the two components. Several contributions will lead to the separation, evolution 
and acceleration of half-solitons. First, their dissociation is triggered by a negative value 
of 02 as discussed previously and can be emphasized by the rolling effect imposed by 
the GL pump. Second, after the unbinding, half-solitons start to feel the longitudinal 
effective field. Let us now discuss the impact on their evolution. 

The normalized condensate wavefunction carrying a dark soliton in both 
components reads 

- _ ' V^+ \ _ [^ ( tanh [/+ {x - d/2)] e^<^" , 
^ " ' V^_ ) " V 2 I tanh [/_ (x + d/2)] ' ^^' 
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Figure 3. (a) Normalized slice of a GL beam together with its phase (for the case 
I = +1) as seen by a pure ID system, (b) Scheme of the separation of a Gaussian 
beam into two GL beams subsequently shined on a wire shaped microcavity. (c)-(f) 
Impact of the traverse displacement of the spots, here a2 = and {1+J-) = (—1, —1). 
(c)-(d) No transverse displacement leads to static dark HSs in each component, (e)-(f) 
A small y shift leads to a uniform propagation of the HS. We see that the direction and 
speed of propagation depends on the direction of the shift sign and on the proximity 
of the spot center from a boundary respectively. In (e) the a+ (ct_) beam is shifted 
by Ay_|- = +0.25 fim (Ay_ = —0.50 /im). The blue lines in (c) and (e) stand for the 
wire's boundaries while the white dashed lines show the position of the spot center. In 
(d) and (f) the white regions are absent of particles due to the finite size of the spot. 



assuming a variable spatial separation d and 0:2 = 0. 0o is a constant relative 
phase between the two components and (Z+, 1-) can take independently the values 
±1 determining the sign of the vr phase shift through each soliton. The pseudospin 
components are straightforwardly calculated using the Eqs.(6) giving 
n , , , , r, / d^ 
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In the case where d = 0, it is seen that the orientation of 5* is homogeneous and fixed 
by the total relative phase A^ = (/)osign(/-|_//_). On the other hand for d 7^ 0, while 
the homogeneous pseudospin texture remains constant far from the solitons cores, it is 
of course modified between them. Comparing e.g. S^^|(±oo) and 5'||(0), it is easily seen 
that the direction of the in-plane pseudospin is always opposite between the solitons 
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and away from them. This is a crucial point since depending on the orientation 
of the effective magnetic field ^lt, the half-soliton separation will occur in opposite 
directions. For example with the set (/+,/_, 0o) = (+1)+1)0), A6 = and therefore 
S^^l(±oo) = +SxUx while S\\{—d/2 — )■ d/2) = —S^Ux- A field Qlt = ^xUx is parallel 



to 5*11 far from the solitons bringing negative contribution in to Emag [see Eq. (IT2l) ]. the 
latter is therefore maximized increasing the spacing between the solitons, where Qlt 
and S\i are antiparallel, the solitons are consequently accelerated in opposite directions 
(because their mass is negative). In this configuration, the integer soliton is unstable 
against the field, the slightest symmetry breaking (e.g. some noise) between the a+ 
and cr_ components result in its decay into half-solitons. The situation is opposed 
for (/+,/_, 0o) = (+l,+l,7r). We finally note that AO = {p + 1/2)tt {p G Z) implies 
S\\{x)-Qlt=0 giving Emag{d) = Eq: An integer soliton or a pair of half-solitons carry the 
same magnetic energy whatever d. However, as soon as one of the half-solitons starts to 
move, the S\i (x) ■ Qlt=0 is no longer verified, and the half-solitons become accelerated 
by the field. Considering separately each half-solitons, their charge can be expressed 
as g = go(l ~ f ^/c^) with go = sign(/) cos(0o) and the divergent monopole-like texture 
appears e.g. for (/, 0o) = (+1, 0). We show in the Fig H] several types of soliton textures 
depending on the value of A9. 



Ae = 
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Figure 4. Pseudospin textures of four separated half-soliton pairs depending on the 
total relative phase A0. 



To underline the contributions to the solitons dissociation, we have numerically 
implemented a realistic configuration using GL beams in each component with winding 
numbers (/+,/_) = (+1,-1-1) and 0o = 0. We show in Figj5] the HS propagation 
separating the different contributions in a situation where there is no rolling effect. 
In the panel (a), we show the 02 impact on the separation, with '^x = 0: we observe the 
linear trajectories (similarly to Figj2]^c)) after the HS are released from the pulsed pump 
spot, no acceleration is observed. In the panel (b) only Qlt = ^xUx is present (^2 = 0), 
trajectories become parabolic up to the limiting speed, which is a clear signature for the 
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constant monopole acceleration (non-relativistic limit). In the panel (c), we show the 
combined effect of both 02 and Qlt together with the traces of the solitons trajectories 
(dashed lines) from the panels (a) and (b), in that case the acceleration is emphasized 
by the initial repulsion. The rolling effect can either assist the separation or block it, 
depending on both the winding number imposed and the lateral shift direction. 




(c) 



50-50 




x{Mm) 



Figure 5. Contributions to tlic HS motion, (a) a2 = —0.2ai and hilx = 0, (b) 
a2 = Oai and hfl^ = 100 /zeV and (c) a2 = — 0.2ai and hil^; = 100 ^eV. The dashed 
lines stand for the trajectories in the configuration of the panel (a) (white) and (b) 
(black) to be compared to the red line. 



We have seen that the splitting of a linearly polarized beam into its circularly 
polarized components allows an independent selection of the winding number of the 
half-solitons in each component and to fix their initial separation. However, even this 
step can be avoided, in order to further simplify the experimental configuration. Indeed, 
making the entire beam scatter on a hologram will impose the same winding number 
/+ = /_ to both condensate components, and the two components will this time be 
perfectly superimposed on the sample. In that case, the separation will not occur, 
unless the symmetry between the a+ and o"_ components is broken by some means. 
Although noise will always be present in the system and allow the separation to occur, 
it might induce (in some realizations) a symmetry breaking on time scales larger than 
the polariton lifetime. Moreover, the 0"+ and o"_ solitons will be separated randomly 
(e.g. cr_ going to left and cr+ going to the right in one realization and the opposite 
behavior in the next one). These two points, that we have checked numerically (not 
shown), are clearly harmful to the reproductibility of the effect, and therefore, to the 
creation of a spin current. One should also bear in mind that a beam prepared to be 
linearly polarized might carry a small ellipticity. The ensuing density imbalance between 
the two condensate components will lead to the formation of HSs with slightly different 
healing lengths and therefore different effective masses. The rolling effect (which will 
occur in the same direction for both HSs) or possibly the wedge naturally present in 
microcavities (or the small gradient of the wire width), will induce an effective mass- 
dependent motion providing the separation. We have simulated this configuration with 
no initial separation, (/+, /_) = (+1, +1) and a 1% ellipticity of the input beam. We show 
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in FigJ7]the separation obtained from the rolling effect shifting the beam by Ay = +0.5 
fim [panel(a)] and from a wedge in the sample producing linear ramp potential of slope 
10 /ieV//im [panel(b)]. In that latter case the HSs are accelerated by the force they 
undergo not from the side of the magnetic field, but directly from the potential {hQ^ = 
here). 




x(Mm) 



x(|jm) 



Figure 6. Separation of perfectly overlapping HSs thanks (a) to the rolling effect 
induced by a Ay = 0.5 /um shift and (b) to a wedge producing a potential ramp of 
slope 10 yueV/yLtm. For both configurations a2 = — 0.2ai, h^x = and 1+ = I- = +1. 

In conclusion, polaritons condensates in quasi-lD microwires constitute an ideal 
system for the creation of magnetic currents based on half-solitons thanks to their 
easy nucleation, natural separation and acceleration. The half-light component in the 
polariton wavefunction allows the particles to travel at high velocities and therefore the 
analogues of magnetic monopoles that are half solitons can travel at the speeds close 
to the speed of light. They constitute extremely promising entities for the fabrication 
of a new class of high speed spin-optronic devices. While the propagation of these 
magnetic charges in wires is technologically attractive, two dimensional systems (planar 
microcavities) allow even more freedom and fundamental richness encompassed in the 
angular momentum of vortices as we will discuss in the following section. 



4. Half-vortices separation 

Integer vortices in "scalar" polariton condensates (that is, circulary polarized, without 
significant coupling between the two components) have already been demonstrated 
experimentally in the optical parametric oscillator (OPO) configuration |25ll23]. These 
vortices have been shown to be relatively persistent, remaining in the ground state of the 
OPO for long times, much longer than polariton lifetime. Another configuration, where 
the integer and half-integer vortices occur due to the persistent flows in the polariton 
condensate under non-resonant pumping, has also been studied experimentally [261 122] ■ 
In this work we are not particularly interested in demonstrating the self- sustained 
coherence of polariton condensate maintained by the OPO. We propose a configuration 
which seems to be the simplest from the point of view of experimental realization, 
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in order to observe controlled separation of half-vortices. Using the Gauss-Laguerre 
beam we create an integer vortex, which then evolves freely under the effect of spin- 
dependent polariton-polariton interactions, constant effective in-plane magnetic field, 
and fc-dependent TE-TM splitting. We demonstrate numerically the separation of a 
vortex (/_!_,/_) = (+1,+1), with both components rotating in the same direction, in 
a weak static in-plane field of 10 /xeV. Such integer vortex is likely to form in case of 
non-resonant pumping, due to persistent flows in the polariton condensate in a disorder 
landscape. The separation into half-vortices occurs for 0o = ^r, while for 0o = the two 
vortices remain coupled with each other. Indeed, the half-vortices of opposite charges 
form a dipole, and the interaction between these charges induced by the field increases 
with the applied field. This is why, depending on the initial winding numbers and on 
the pseudospin texture defined by 0o, one can sometimes observe the stabilization of 
the dipole length at an equilibrium value. Figure 7 shows the initial pseudospin texture 
(a) and the comparison between the situations without in-plane magnetic field (b) and 
with magnetic field (c). In the case (b) the small separation is due to the a2-induced 
repulsion. In the panel (c) the half-vortices are dissociated and accelerated in opposite 
directions by the constant in-plane magnetic field, although its value is relatively small. 
The parameters of the Gauss-Laguerre beam have to be chosen carefully, in order 
to minimize the non-desired effects such as the dynamic formation of solitons on the 
outer horizon of the density profile. At the same time, the size of the minimum in the 
center should be in agreement with the expected healing length. In a word, one should 
be as close to the perfect initial condition of a flat infinite condensate with a vortex in 
the center, as possible. 

™ — — — ///^^' — "' — i^-. ^^C wwww-^ \ t / 

_, — ____., — __-^ z/^^^--^— — !rv\\ "^"^^"^^ w /-"^v \ / 
^ . _^ :::: ,W»\\\ \ )iC\ \ /^^::i.(P'/ 

__^ C^^-^\ t m \ \ t I iC^^ \ f /^^^^^^^ 

_< _. .::_^v\ \ \ \ \ \ \ \ ] c^\ \ f /^^^^^^^ 

on — ' -- -^. -"-^ \ \ \ \ M f ; _»^ \ t / / ^ ^^^y^ y 

-ZU , , __^-.\\ \ \ \ M f , .'^\ \ f / y yyyyy/' 

-20^.0 , 20 -20^,0 .20 -20^ .0 . 20 
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Figure 7. a) Initial pseudospin configuration for an integer vortex (+1,+1); b) 
Dissociation into two HVs under tlie effect of a2-iiiduced repulsion only; c) Dissociation 
and acceleration of the two HVs under the effect of magnetic field 17^;. Panels (b) and 
(c) correspond to 50 ps after the pulse. 

However, peculiar undesired effects can arise if a half-vortex becomes trapped in 
the spatial density inhomogeneity created by the rotation of the pseudospin of the other 
half-vortex, which is not cylindrically symmetric: in some directions the pseudospin is 
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initially aligned with the field and does not rotate, while in the other directions the 
pseudospin is not aligned with the field and is bound to gain a nonzero Z-projection, 
which can be a barrier or a trap for the other half-vortex. 

When the half-vortices are sufficiently far from each other (that is, when the 
distance between the vortices d :» ^), each of them can be considered separately, as 
a half-vortex composed of a vortex in one component and a homogeneous background 
in the other component. Each of them is therefore subject to a force acting from the 
effective magnetic field, and we can expect them to accelerate freely. The pseudospin 
texture of a half- integer vortex can be either convergent or divergent, looking similar to 
that of a point charge, as it is the case in ID, but there are two important differences 
between the ID and the 2D case. Indeed, in ID the constant pseudospin field far from 
the vortex core is exactly the solution of the Maxwell's equation V ■ S = S{x), whereas 
in 2D the pseudospin texture does not depend on the distance from the core (as well as 
in ID), while the solution of the Maxwell's equation in this case should be decaying: the 
field of a charge depends on the distance from it. The second important difference resides 
in the presence of the relative phase (pQ, which influences the more complicated textures 
of the half-vortex and its propagation direction. Various cases and the corresponding 
textures have been considered in Ref. [3]. Moreover, in the 2D case there is a long- 
distance interaction between the vortices lying in the same component (which is not 
the case for solitons in ID), and even between vortices in different components, when 
a magnetic field is applied. The field creates a transfer of particles between the two 
components, and since the particles have nonzero propagation velocities even far from 

— * 

the vortex core, this creates a long-distance interaction force, absent in ID. For the k- 
dependent TE-TM splitting in planar microcavities, this interaction has been considered 
in Ref. [28]. 






*^J^A 







Figure 8. Pseudospin textures of (a) phase vortex (Z4., Z_, (/jq) = (+l!+li''') and 
(b) a polarization vortex (/+,?_, 0o) = (1,-1, tt), dissociated in their associated half- 
vortices. 



In the ideal case, the condensate wavefunction containing a half-vortex in both 
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components separated by d reads 




n_e 



The radial functions take the approximated form pj n±{r±,(l)±) = rj^/{rj^ + 2), where 
r± = -\/(x~^~rf72p~+T^~TW2p and 0± = arctan[(?/^(i/2)/(x=F'^/2)]. The components 
of S are easily found using the Eqs.(6). To give the complete picture, we show in 
FigJS] the pseudospin textures for two half-vortex pairs using the sets of parameters 
(/+,/_, 0o) = (+l,+l,7r) and (/+,/_, 0o) = (1,-1, tt) which constitute a phase and a 
polarization vortex respectively when the half- vortices overlap. 

In conclusion, half-vortices embody the 2D candidates for the monopole analogy 
possessing a divergent pseudospin field (depending on 0o)- They accelerate under the 
action of the effective field and are characterized as well by a 0o-dependent propagation 
direction [3] at fixed orientation of the effective field. They constitute the building 
blocks for potential 2D-magnetic circuits. 

5. Conclusions 

We have shown that the half-integer topological defects in spinor polariton condensates 
not only possess a divergent or convergent pseudospin texture similar to the textures of 
point charges in ID and 2D, but also behave as magnetic charges in presence of effective 
magnetic fields, accelerating along them. We have proposed a simple configuration for 
experimental creation, dissociation, and acceleration of half-integer topological defects. 
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